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We theoretically study the optical properties of a Fermi- 
Dirac gas in the presence of a superfluid state. We calcu- 
late the leading quantum-statistical corrections to the stan- 
dard column density result of the electric susceptibility. We 
also consider the Bragg diffraction of atoms by means of 
light-stimulated transitions of photons between two intersect- 
ing laser beams. Bardeen-Cooper-Schrieffer pairing between 
atoms in different internal levels magnifies incoherent scatter- 
ing processes. The absorption linewidth of a Fermi-Dirac gas 
is broadened and shifted. Bardeen-Cooper-Schrieffer pairing 
introduces a collisional local-field shift that may dramatically 
dominate the Lorentz-Lorenz shift. For the case of the Bragg 
spectroscopy the static structure function may be significantly 
increased due to superfluidity in the nearforward scattering. 
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I. INTRODUCTION 

As a result of dramatic progress in cooling and trap- 
ping of alkali- metal atomic gases the quantum statistics 
of atoms has observable consequences. Perhaps, the most 
striking effect is the Bose-Einstein (BE) condensation 
of bosonic atoms jj]-|3|- Nevertheless, Fermi-Dirac (FD) 
gases are also expected to exhibit a rich and complex be- 
havior. One especially fascinating property of FD gases is 
that with effectively attractive interaction between differ- 
ent particles the ground state of the system may become 
unstable with respect to the formation of bound pairs of 
quasiparticles or Cooper pairs iQjEJ. This effect is analo- 
gous to the Bardeen-Cooper-Schrieffer (BCS) transition 
in superconductors. FD gases have been an active sub- 
ject of research already quite some time Jfj|-p2[. 

We recently showed 17 that the coherent quasiparti- 
cle pairing between atoms in different internal levels may 
enhance the optical interactions. In particular, the reso- 
nance line of a FD gas (the extinction of light from coher- 
ent laser beam) is broadened and shifted as a result of the 
BCS pairing. In this paper we present a more detailed 
study of the propagation of light in a dilute FD gas in the 
presence of a superfluid state. We also investigate the 
signatures of the BCS pairing in light-stimulated tran- 
sitions of photons between two intersecting nonparallel 
laser beams. We demonstrate the dramatic dependence 
of the Bragg diffraction rate on the BCS order param- 
eter. The Bragg diffraction of atoms has already been 
experimentally used as a method to split a BE conden- 
sate f25f| and to perform spectroscopic measurements on 
the condensates p3,55|. 



In this paper we theoretically study the optical re- 
sponse of a superfluid state in a zero-temperature FD gas. 
First, we consider the propagation of low- intensity light 
and calculate the leading low-density correction to the 
standard century-old column density result of the elec- 
tric susceptibility, also known as the Clausius-Mossotti 
relation po| . This correction is a direct consequence of 
the quantum-statistical position correlations between dif- 
ferent atoms that modify the optical interactions at small 
interatomic separations. 

FD statistics exhibits a short-range ordering of atoms 
in the gas. In the absence of a superfluid state the FD 
statistics forces a regular spacing between the atoms in 
the same internal state within the characteristic correla- 
tion length £f j ~ 1/kp. Here kp denotes the Fermi wave 
number. Fermionic atoms repel each other and short 
interatomic separations are suppressed by the Pauli ex- 
clusion principle. The pair correlation function between 
two atoms displays antibunching. 

The resonant dipole-dipole interactions between differ- 
ent atoms, which behave as 1/V 3 , are dominant at small 
interatomic separations. The Fermi repulsion suppresses 
these optical interactions and therefore the incoherent 
scattering of light in the atomic sample is reduced. As 
a result a FD gas exhibits a dramatic narrowing of the 
absorption line for the coherently scattered light p5| , ^6[ . 

In the s-wave BCS transition the particles near the 
Fermi surface having opposite momenta and different in- 
ternal quantum numbers tend to appear in pairs. This 
leads, e.g., to a finite energy gap A in the excitation 
spectrum of the system and to a nonvanishing expecta- 
tion value of the anomalous correlation function for the 
matter-field annihilation operators (ip-\ (r)^|(r)). Here 
the two internal states are referred to as | and |. 

The BCS s-wave pairing also introduces short-range 
ordering of atoms within the characteristic correlation 
length ~ ep/(Akp) between atoms in different inter- 
nal states 0]. Here ep denotes the Fermi energy. Due 
to the macroscopic two-particle coherence the atoms in 
different internal levels attract each other and the BCS 
pairing enhances small interatomic separations. The pair 
correlation function between two atoms with different 
quantum numbers displays bunching. As a result the 
optical interactions and the incoherent scattering of light 
in the atomic gas are enhanced. This broadens the ab- 
sorption line for the coherently scattered light. 

As a second topic we study the Bragg spectroscopy of 
a FD gas with BCS pairing. In the Bragg spectroscopy 
an optical potential couples to the local number density 
of atoms and connects two external atomic states of dif- 
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ferent momenta. The scattering rate of atoms is propor- 
tional to the dynamical structure function of the atomic 
gas [p pi|j25[ |. The structure function has been exten- 
sively studied for the case of incoherently scattered light 
in optically thin BE |,||-|| and FD atomic 
gases. 

The structure function contains distinct quantum- 
statistical features in the case of fermionic and bosonic 
atomic gases 0. For a FD gas it demonstrates the 
Fermi inhibition of the spontaneous scattering of pho- 
tons P, pl)||T4| ]. The structure function may also be used 
as a method of determining the relative phase between 
two BE condensates [|8|j30j. 

In this paper we show that a superfluid state may sig- 
nificantly increase the value of the static structure func- 
tion. This is because due to the BCS pairing atoms and 
holes near the Fermi surface are mixed; with a given recoil 
momentum there exist more unoccupied states to which 
atoms can scatter. The effect of superfluidity is stronger 
for nearforward scattering. 

In Sec. H we introduce the Hamiltonian density for 
ground-state atoms in the absence of the driving light 
field. The analysis of the quasiparticles follows the stan- 
dard BCS theory Q. We emphasize the effect of the 
quantum statistics on the pair correlation function. In 
Sec. jTj we study the propagation of light in a FD atomic 
gas. The interaction between light and matter is dis- 
cussed in general terms in Sec. [II A . The atomic polar- 
ization is solved for the low-density atomic gas, complete 
with the dependence on the atomic level scheme and on 
various light polarizations in Sec. 1IIB. The electric sus- 
ceptibility representing the damping and the phase ve- 
locity of the light beam is obtained in Sec. 



QIC. The 



effects of the quantum statistics and the s-wave interac- 
tions manifest themselves in two distinguishable param- 
eters. A few rema rks about the light propagation are 
made in Sec. Ill D . In Sec. we consider the possibili- 
ties to observe the BCS pairing via Bragg spectroscopy. 
The Bragg diffraction of atoms in the la ser fie ld probes 
the structure function of the gas. In Sec. [IV B we calcu- 
late the structure function for a superfluid state. A few 
concluding remarks are made in Sec. [v| 



II. GROUND-STATE ATOMS 

A. Hamiltonian 

We assume a FD gas occupying two different internal 
sublevels \g, |) and \g, |) of the same atom with electron- 
ically excited levels |e, v). In the absence of the driving 
light field, atoms in the electronic ground state are de- 
scribed in second quantization by the Hamiltonian den- 
sityW 9 g|: 



Tin 



(1) 



where tp gu (rt) denotes the atom- field annihilation op- 
erator for level \g, v) in the Heisenberg picture, [i gv is 
the corresponding chemical potential, and H-! v m stands 
for the center-of-mass (cm.) Hamiltonian. We have 
approximated the finite-range interparticle potential by 
a contact interaction with the strength given by u g = 
\ixaghjm. Here a g and m denote the s-wave scattering 
length and the mass of the atom. The atoms in differ- 
ent internal states can interact via s-wave scattering. On 
the other hand, due to the Pauli exclusion principle there 
only is a very weak p-wave scattering between two atoms 
in the same level, which is ignored in Eq. (111). 



B. BCS pairing 



Before the light is switched on, the system is described 



by the Hamiltonian density TL = TL g [Eq. 



The as- 



sumption that the driving light only weakly disturbs the 
system allows us to evaluate the electric susceptibility 
by using the ground-state atom correlations determined 
by Tig, even in the presence of the driving light. We 
assume a homogeneous atomic sample and introduce a 
plane- wave basis for the field operators: 



Vv( r ) = -Jy^2 b ^ el 



(2) 



In the Hamiltonian (jl]) we introduce the standard canoni- 
cal transformation to the Bogoliubov quasiparticles plS 



/3-k = "kfr-kT + «kfr ki , 



(3a) 
(3b) 



where u k and w k are real, depend only only on |k|, and 
satisfy u 2 , + v k = 1. The requirement that linearized 
mean field fluctuations of TL g in Eq. ([!]) be diagonal in the 
quasiparticle representation sets an additional constraint 
and we obtain 



(4) 



where E k = yj A 2 + ££, £ k = e k - p + hu g {pi + p{)/2, 
and the energy gap 



h u 



A = -^$>«k(l 



(5) 



In equilibrium, the quasiparticle occupation numbers 
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(/3 k /3 k ) satisfy FD statistics 
^ e E k /k B T + -^-l The dispersion re- 
lation for free particles is given by e k = h 2 k 2 /(2m) and 
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the average of the chemical potentials is /2 = (fx-f +//J/2. 
For simplicity, we assume /i-f = /ij.. For the gap param- 
eter at T = we use the weak-coupling approximation 
A~ 1.76fc B T c |i, where 



k B T c 



8£f 7 _2 
e ' exp 

7T 



(6) 



Here the Fermi wave number kp = (6tt 2 p) 1 ^ 3 is defined in 
terms of the atom density p, T c denotes the critical tem- 
perature of the BCS phase transition, and 7 ~ 0.5772. 

In the superfluid phase transition the atoms in the 
different internal states f and J, form quasiparticle 
pairs resulting in a nonvanishing anomalous correlation 
(i/'|(ri)'0x( r 2))- The effect of this macroscopic two- 
particle coherence on the atomic position correlations is 
particularly transparent in the case of the pair correlation 
function: 

p 2 (ri^,r 2 crr) = (^ iy (r 1 )^| CT (r 2 )^ gT (r2)'0 ff7) (r 1 )) . (7) 

In the ground state of TL g [Eq. (Q)], determined by the 
vacuum of the Bogoliubov quasiparticles [Eq. (||)], the 
pair correlation function reads (for v 5^ a) 

p 2 {r 1 w,T 2 (J(T) = p v p a + \(ipgu{ri)tp ga (r 2 ))\ 2 , (8a) 

p 2 (v 1 w,v 2 vv) =pl- |(V^(ri)T/V(r 2 ))| 2 . (8b) 

The pair correlation function represents the joint proba- 
bility distribution for the positions of two atoms. We note 
from Eq. (8b) that the fcrmionic atoms in the same inter- 
nal state repel each other analogously to antibunching of 
photons and short interatomic separations are inhibited. 
For an ideal homogeneous FD gas, in the absence of a 
superfluid state, we can analytically evaluate Eq. (|Sb|). 
We obtain 111 



p 2 (r\vv, vv) = p 2 



sin kpr 
kpr 



cos kpr 



(9) 

where r = \v\ — v 2 \. In Fig. [l] wc have shown the pair 
correlation function (^). The FD repulsion between dif- 
ferent atoms results in suppressed dipole-dipole interac- 
tions and in the Fermi inhibition of incoherently scattered 

light n@. 



p 2 /p 2 




FIG. 1. The pair correlation function for an ideal homo- 
geneous Fermi-Dirac gas at zero temperature. The two atoms 
have the same internal quantum numbers. The pair correla- 
tion function represents the joint probability distribution for 
the position of the second atom given that the first atom sits 
at the origin. Two fermionic atoms repel each other and small 
interatomic separations are suppressed. 

As a result of the BCS pairing atoms in different in- 
ternal states attract each other analogously to bunching 
of photons according to Eq. (|8a|). Therefore short inter- 
atomic separations arc enhanced. In the next section we 
find that this may lead to enhanced optical interactions 
and incoherent scattering of light. 



III. LIGHT PROPAGATION 

A. Basic relations 

In this section we introduce terms in the Hamilto- 
nian that result from the electromagnetic fields. Wc 
briefly recapitulate and generalize our previous quantum 
field-theoretical analysis of the light-matter interactions 
0'H 15 0J- As a basic assumption, atoms are repre- 
sented as point dipoles and the radiated field in a m ediu m 
has the familiar expression of the dipolar field [Eq. (|15f 



1. Electromagnetic field 

We consider the nonrelativistic Hamiltonian formal- 
ism of electrodynamics. It is advantagous to study the 
propagation of light by introducing the dipole approxi- 
mation for atoms and the corresponding Hamiltonian in 
the length gauge obtained in the Power-Zienau-Woolley 
transformation [f^-f^s] . 

The electromagnetic field introduces additional terms 
in the system Hamiltonian. The Hamiltonian for the free 
electromagnetic field-energy is 



d rTLp 



^ huj q a\a q . 



(10) 



Here u> q and a q denote the mode frequency and the pho- 
ton annihilation operator. The mode index q incorpo- 
rates both the wave vector q and the transverse polariza- 
tion e q . In the length gauge the basic dynamical degree 
of freedom for the light field is the electric displacement 
D(r) which interacts with the atomic polarization P(r), 



H D = --P(r) -D(r). 



(11) 



In the dipole approximation the polarization is given in 
terms of the density of atomic dipoles P(r) = J2i di<5(rj — 
r). Here and denote the dipole operator and the 
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cm. position operator for the ith atom. In second quan- 
tization the positive frequency component of the polar- 
ization is given by 



P+(r) - 



£ d s!/e7) V|,(r)Ve^(r) = P U r ) > ( 12 ) 



I'.i) 



Vyf) 



p ^( r ) = dg^ ev ^l„(r)ilJe V (r) , (13) 

where d gveri stands for the dipole matrix element for the 
transition | e, 77) — > \g, v) 



V^2e a (eij; lg\l<r;gv) 



(14) 



Here T> denotes the reduced dipole matrix element, which 
is chosen to be real, and d er]gl , — A* guerj . The summation 
in Eq. ([lj]) runs over the unit circular polarization vectors 
(j = ±1,0 weighted by the Clebsch-Gordan coefficients 
of the corresponding optical transitions. Here light fields 
with the polarizations e± and eo drive the transitions 
\g, v) — > \e,v ± 1) and \g, v) — > |e, v), respectively. 

For a weak external magnetic field the nuclear spin /, 
the electron spin S, and the orbital angular momentum 
L may be coupled. In that case F and mp, defined by 
F = I+S+L, are good quantum numbers. With stronger 
magnetic fields the nuclear spin decouples and the optical 
transitions are between different sublevels of J, where 
J = S + L. In the latter case we may assume that the 
atoms occupy a single manifold of to/. In the following 
calculation we consider atomic sublevels determined by 

TO p . 

The positive frequency component of the electric field 
E + may be expressed in terms of the positive frequency 
component of the driving electric displacement D^., with 
the wave number fc, and the dipole radiation field |3l|: 



e E+(r) = D+(r) + J d 3 r> G(r - r') P+(r') , 



(15a) 



Gij(r) 



' d__d_ 

dri drj 

We also have the familiar relation 



ikr 



e" 



D(r) =e E(r)+P(r). 



M(r) ■ (15b) 



(16) 



where n = r/r. The volume integral over 1 / r 3 in Eq. 
is not absolutely convergent in the neighborhood of the 
origin. The expression ( |l7| ) should be understood in such 
a way that the integral of the term inside the curly brack- 
ets over an infinitesimal volu me e nclosing the origin van- 
ishes J3lj . We note that Eq. (|l5a|), with the correct delta 
function contribution from Eq. (|17|), yields V • D(r) = 
for neutral atoms justifying the use of the electric dis- 
placement, instead of the electric field, in Eq. (fil"|). 

The nonrelativistic propagator ( [l7|) involves an explicit 
high-frequency cutoff pl| . In situations where integral 
expressions containing the propagator G(r) are not ab- 
solutely convergent the integrals are defined in such a 
way that the integration over spherical angles should be 
performed first |pl[| . 

In the length gauge the Hamiltonian also contains the 
polarization self-energy term: 



2e 



P(r) • P(r) 



(18) 



However, this is proportional to the overlap of the atomic 
dipoles and, in the limit of low light-intensity, all contact 
interaction terms between different dipole atoms were 
shown in Ref. to be inconsequential for light-matter 
dynamics. Although the result is valid for an arbitrary 
number of atoms, the cancelation of the contact interac- 
tions is especially transparent in the case of a superradi- 
ant decay of only two point dipoles [B2J . 



2. S-wave interactions 

Two cold fermionic atoms in different internal states 
interact by means of s-wave scattering. This introduces 
interaction terms between ground-state atoms in Eq. (jl]) . 
In the presence of driving light we also have s-wave scat- 
tering between different electronically excited levels and 
between electronically excited and ground-state atoms. 

We write the contribution to the Hamiltonian density 
that consists of the excited-level operators as 



(19) 



Equation (|l5|) is the integral representation for the 
Maxwell's wave equation of the electric field in the pres- 
ence of dipole atoms and the monochromatic dipole radi- 
ation kernel G(r) coincides with the corresponding clas- 
sical expression [p6[ . The explicit form of the radiated 
field from a dipole with the amplitude p reads: 



G(r)p 
+ [3n(n • p 



fc 3 f e 



P][ 



1 



Akr 



{krf {krf 



p(5(r) 



(17) 



Here u e 



4:Trha P 



-/to describes the two-body 



s-wave scattering between the atoms. For simplicity, the 
frequency of the optical transition ujq is assumed to be 
independent of the atomic sublevel. For typical values 
of the optical lincwidth the cm. motion for the excited 
atoms may be omitted ||. In this paper we consider a 
situation where the intensity of the driving light field is 
low and therefore the density of the excited atoms is low. 
Hence, to leading order in the limit of low light-intensity 
the second term in Eq. (|19|) makes no contribution to the 



4 



optical response, and we do not address its explicit form 
in more detail. 

We assume that to leading order all remaining interac- 
tions between the ground-state and excited-state atoms, 
which cannot be accounted for when the atoms are mod- 
eled as point dipoles, are governed by the following in- 
teractions p5| : 

Kge = ft ^2 Ugunearlpl^t^e^gT • (20) 

If we assume the conservation of the angular momentum 
of the colliding particles, the two-body interaction can 
be written in the following form p6[ : 

fa+fb f 

n ge = h J2 E u F , mF FmF Fmp , (21) 

F=\f a -f„\m F =-F 

where UF,m F = ^ft<iF,m F /tti and 

Fmp = ^2 { Fm F\fah\fa,m a ]f b ,m, b )ll)f ama 7p fbmb . 
m a mh 

(22) 

Here (Fm F ; f a f b \f a ,m a ] f b ,m b ) is the Clebsch-Gordan 
coefficient for forming the state \Fm F ) from two colliding 
particles in the angular momentum states \f a m a ) and 
\hm b ) . 

As an example we consider the / = 1/2 —> 3/2 transi- 
tion. By expanding Eq. (^Tj) we obtain H ge — Hfe+Hge, 
where 7i^ e does not mix the population between the dif- 
ferent sublevels. The explicit form reads 

H ge = h H u gv^\A\ v ^aA™ ■ (23) 

For the / = 1/2 — > 3/2 transition the summation runs 
over v = ±1/2 and cr — ±1/2, ±3/2. The interaction 
strengths u gvfia are determined in terms of the interac- 
tion strenghts u Ft m F for the scattering channels \Fm F ): 

U g±i,e±| = U 2,±2, 

U g± i eT i = -(ui.o + U2,0), 
3 1 

1 3 

U S ±i,e±i = 4«l,±i + -U2,±l. 

On the other hand, Hg e consists of spin-exhange collision 
terms: 

+ —("2,1 - Ul,l) ^|^_i^ ff i^ei 

+ —("2,-1 - U lt _l) ^e-|^| + RC - ■ ( 24 ) 



Similar nonlinear wave-mixing terms can have interesting 
effects on the dynamics of spinor BE condensates |37| . 

If the s-wave scattering amplitude u F , mF is indepen- 
dent of F we obtain Hg e = 0. Furthermore, if the scat- 
tering is also independent of m F , u gv ^ ea in Eq. ( p3| ) docs 
not depend on sublevels v and a . For simplicity, in the 
following we assume that Hf e =0. 

B. Atomic polarization 

The dipole radiation [Eq. (fL5|) l describes the scattered 
light in a medium. In this section we study the optical re- 
sponse of an atomic gas by solving nonperturbatively the 
polarization of the matter-field for a low-density atomic 
gas in the limit of low light-intensity. The main items 
are the steady-state equation of the atomic polarization 
[Eq. (p9|)] and its solution by means of the low-density 
decorrelation approximation [Eq. (|32|)]. 

1. Equation of motion 

We consider the optical response of the atomic po- 
larization in the limit of low light-intensity. We derive 
the corresponding Heisenberg equations of motion for 
the matter-field operators. Alternative approaches have 
used, e.g., the Schwinger-Keldysh techniques [Q. We 
assume that the spin-exchanging collisions between the 
ground-state and the excited-state atoms vanish indicat- 
ing Hge = ■ Then the Heisenberg equations of motion 
for atomic field operators are obtained from the Hamilto- 
nian densities 7L [Eq. (§], H D [Eq. @], H e [Eq. @], 
andtt£ [Eq. (§|)]. 

In the limit of low light-intensity we derived from the 
Heisenberg equations of motion a hierarchy of equations 
for correlation functions involving atomic polarization 
and atom density (3^^]. In the case of the present 
system we may proceed similarly. As far as the opti- 
cal response is concerned it is again assumed that we can 
concentrate on the dynamics of internal degrees of free- 
dom for the atoms and the light. Hence, in the equation 
of motion for the atomic polarization the kinetic energy 
of the atoms is neglected. Nevertheless, the quantum 
statistics of the different cm. motional states is still fully 
included. 

The vacuum electromagnetic fields are eliminated by 
performing the field-theory version of the Born and 
Markov approximations (3^]. We then obtain the equa- 
tion of motion for the polarization operator component 

P£,(r). 

j t P+ = (iS - 7 )P+ + • D+ 

+ inf dV P^G'(r-r')4,P + (r')Vv 
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in J d 3 r' ?Z ■ G'(r - r')iP+(r')i 

i^L^PeuKl ■ + iu g (l - 5 v>a ) ^g^l^ga^er, 



-\- illgi/^ ecr ^ gv V^ecr V*e<7 ^^go , ,e7]^Pgi/^Pga'^9<^'^GT} 

+ i{u 

er * 



(25) 



Here we use the notational convention that the repeated 
indices er, r, and £ indicate the summation over the corre- 
sponding sublevels. In Eq. ( |2q ) we have shown explicitly 
only the nonlocal position dependence. The atom-light 
detuning is denoted by S and the Kronecker delta func- 
tion by 5 V ^. The spontaneous linewidth 7 is given by 
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V 2 k 3 
67r?ieo 



(26) 



We have also defined k = V 2 /(?ieo) in terms of the re- 
duced dipole matrix element T>. In Eq. ( |25| ) we intro- 
duced the following projection operator 

VI] d.g Ue7 ](l e £g T 



field amplitude. Therefore, we take the low-intensity 
limit by retaining only those products of operators in 
Eq. J25| ) that involve at most one excited state field op- 
erator or the electric displacement amplitud e |3l| . Then, 
e.g., the third line and the last line in Eq. (|25|) make no 
contribution to the equation of motion for P~(r) for low 
light-intensity. 

We also note that in the low-intensity limit the pair 
correlation function (^) is determined by the Hamiltonian 
density in the absence of the driving light H g [Eq. (||)]. 
This is because the effect of the light on the ground-state 
field amplitudes involves terms that contain at least one 
excited-state field operator and one light-field amplitude: 

ijj g (r) cx d ge ■ D~(r)^ e (r) . 

Thus, to leading order in the low-intensity limit the effect 
of the driving light on the ground-state atom correlation 
functions vanishes. 

For the expectation value of the polarization we use 
the notation = (P^), with v and 77 denoting the 

atomic sublevel. The steady-state solution of V\ vri in the 
limit of low light-intensity is given by 



£ H^M{<; MW) , (27) p ^ (ri) = ap ^ . D + (ri) + £^p 2(ri<T<r;riI ^ 

+« £ / ^ P ^ ' G '( ri - r 2 )P 2 (ri^; r 2 r£) . (29) 

<tt(, 



to include the dependence of the scattered light on the 
polarizations and on the atomic level structure. 

In the present formalism the atoms are represented 
by ideal point dipoles which may overlap. Obviously, 
real atoms with a hard-core interatomic potential cannot 
overlap. We remove the contact dipole-dipole interac- 
tions between different atoms. This is done by introduc- 
ing in Eq. (|25|) the propagator 



Here a = —V 2 /[heo(8 + ij)] is the polarizability of an 
isolated atom and p v = (tp^ipg^) denotes the ground- 
state atom density in level v. We have also defined 



P 2 (ri^77;r 2 (TT) = (^t / (ri)P+.(r 2 )V> OT (ri)) 



3,(r)=G«(r) 



<%<5(r) 



(28) 



-par] 

Til 



that explicitly cancels the contact interaction of G(r) dis- 
played in Eq. (|l7|). The purpose of this definition is to 
yield a vanishing integral for G'(r) over an infinitesimal 
volume enclosing the origin. As shown for the case of 
low-intensity light in Ref. we can make the substi- 
tution ( p8|) without changing the outcome of the opti- 
cal response, and without any additional assumption of 
the hard-core interatomic potentials, even for ideal point 
dipoles. This is because due to the divergent dipole- 
dipole interactions all correlation functions for atomic 
dipoles vanish whenever two position arguments are the 
same. As a result the contact interaction terms between 
different dipole atoms do not have any effect on the light- 
matter dynamics. The independence of the optical re- 
sponse of the collection of dipole atoms on the substitu- 
tion ( p8| ) is the underlying explanation for the Lorentz- 
Lorenz local-field correction to the electric susceptibility. 

We consider the limit of low intensity of the driving 
light. Obviously, light has to be present in order to pro- 
duce population in the electronically excited levels and an 
excited-state field amplitude is proportional to the light 



5 + i~/ 



Lgaerj (1 ^TI/)^(/] 



(30) 



(31) 



The normally ordered expectation value P 2 (rizyr;; r 2 o"r) 
describes correlations between an atomic dipole at r 2 and 
a ground-state atom at ri. In the integral of Eq. (pjl) 
it represents a process in which an excited atom at r 2 
emits a photon and excites a ground-state atom at rj.. 
The tensor J 7 ^ generates the collisionally-induced level 
shifts. 



2. Low-density approximation 

So far, we have obtained a steady-state solution for 
the atomic polarization ( p9|) tha t acts as a source for the 
secondary radiation in Eq. (15a). Equation (|29] ) involves 
unknown correlation function P 2 . Basically, we could 



continue the derivation and obtain the equations of mo- 
tion for P 2 and for the higher-order correlation functions. 
This would eventually result in an infinite hierarchy of 
equations analogous to the equations in Ref. |5lJ. How- 
ever, even in the case of a simple level structure and in 



G 



the absence of the s-wave interactions the solution for the 
whole system by stochastic simulations is demanding on 
computer time Q . In the studies of the refractive in- 
dex of a quantum degenerate BE gas Morice et al. [ j39|]40| 
considered a density expansion in terms of the number 
of atoms repeatedly exchanging a photon by introducing 
certain approximations to the ground-state atom corre- 
lations. Although the lowest-order density correction to 
the susceptibility of a zero-temperature FD gas may be 
obtained analytically , in the presence of highly non- 
trivial quantum statistical position correlations a rigor- 
ous density expansion is in most cases a very challenging 
task. In this paper we consider low atom densities (in 
terms of p/k 3 ) and approximate Eq. ( p9[ ) by the decou- 
pling that is analogous to the lowest-order correction in 
Ref.il, 



the vector components of Pi for the / = 1/2 — > 3/2 tran- 
sition having the electronically excited sublevels mt = 
±1/2, ±3/2. The pair correlation function in Eq. ( p2| ) is 
nonvanishing only with v = rj. Because we are dealing 
with a linear theory, the electric field and the polarization 
are related by the susceptibility as P + = 6qxE> + ■ We 
consider a situation where FD gas fills the half-infinite 
space z > 0. For simplicity, we assume equal and con- 
stant atom densities for the spin states p-|- = pi = p. To 
simplify further, we assume that scattering length a gvea 
is independent of v and a corresponding to the case that 
the s-wave interactions in Eq. (|2^) are independent of 
the scattering channel IFmp). We write the incoming 
free field as a plane wave 



D F (r) = D F ee lkz 



(33) 



P^r^ ^^^W fo), (32) 

Pa 

where the ground-state pair correlation function p 2 is de- 
fined by Eq. (@). 

The decorrelation approximation ( |32"|) introduces the 
lowest-order correction to the optical response in terms 
the number of microscopic optical interaction processes 
between the atoms by ignoring the repeated scattering 
of a photon between the same atoms [[it)) . As shown in 
Rcf. Jl6| in the absence of a superfiuid state it also cor- 
rectly generates the leading low-density correction. The 
predictions of the expansion by Morice et al. (39) were 
tested for a zero-temperature FD gas in one dimension 
|jl5f . The agreement with the exact solution obtained by 
the numerical simulations was found to be semiquantita- 
tive and in the low-density limit excellent. 

The dependence of the light propagation on the den sity 
fluctuations may now be observed by inserting Eq. ( |32j ) 
into Eq. ((2^). If the emitting atom at r 2 and the ab- 
sorbing atom at i"i have the same internal state, the pair 
correlation function displays repulsion and is determined 
by Eq. (Sb). In the case of different internal states the 
atoms attract each other and the pair correlation func- 
tion is obtained from Eq. (|Sa|). 

It is crucial for the low-density limit that the atom 
operators in Eq. ( p5| ) were arranged to normal order. 
Otherwise, commutators are generated for higher-order 
correlation functions P2 , • ■ ■ that could be of the same 
order in atom density as the terms in the equation for 

Plvrj' 



C. Electric susceptibility 

In the previous section we obtained the steady-state 
solution for the atomic polarization ( p9| ) by means of the 
low-density approximation (p2|). The optical response 
may now be evaluated by eliminating and P2 from 
Eqs. (15a), (p9l), and (|32]). As an example we calculate 



We assume that it is linearly polarized with e parallel to 
d g i e i . By choosing e = e x we have the following repre- 
sentation for the circular polarization vectors in terms of 
the unit Cartesian coordinate vectors 



1 



eo 



(34) 



With the ansatz Pi vv (r) = Peexp(ik'z), for Im(fc') > 
0, we then immediately see that Pi ur) — for v 7^ 77. 
Finally, by using Eq. (g), and by ignoring the effects of 
the surface of the atomic gas fL5| , we obtain a spatially 
constant susceptibility for the sample as 



with 
Si = 



X ~ k 2 



Ca 

/ (fre 

P 



2Cap 



1 - 2Cap/3 + Hi + S 2 



(35) 



—ikz 



G'(r).e[KV sT (r)V si (0))| 2 



■K^(r)^l(0))| 



(36) 
(37) 



Here we have used the obvious relation p2{v\a, y 2 v) = 
P2{y\<J<J^ivv) — pi{v\vv,Y%a<j). The parameter C de- 
notes the value of the Clebsch-Gordan coefficient C — 
|(ei;l(7|10; 5 i)| 2 = 2/3 in the case of the / = 1/2 -> 3/2 
transition. By writing r = kr the propagator in Eq. ( |3q ) 
has the following expression in the spherical coordinates 



e* -G'(r)-e x = 



k A e 



47T 



(1 — sin 2 9 cos 2 



1 



±(3sin 2 6»cos 2 0-l)(— 



(38) 



In Eq. (|3^) Si is solely generated by the quantum- 
statistical position correlations between different atoms. 
The effect of s-wave interactions is encapsulated in £2 • In 
an uncorrelated atomic sample the atomic positions are 
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statistically independent and the pair correlation func- 
tion (0) satisfies p%{rv : r'a) = p v p a resulting in Si = 0. 
For the case of uncorrelated atoms, and in the absence 
of the s-wave scattering, we would obtain Eq. (]3^) with 
Si = S2 = 0. This is the standard column density re- 
sult stating that susceptibility equals polarizability of an 
atom times atom density. Equation ([35]) also involves the 
Lorentz-Lorenz local-field correction in the denominator. 

The quantum-statistical corrections to the column den- 
sity result are introduced by Yi\. It describes the modi- 
fications of the optical interactions between neighboring 
atoms due to the position correlations. The second term 
in Eq. (^6|) represents the quantum-statistical contribu- 
tion to the scattering process in which a photon emitted 
by an atom in internal state v at position r is reabsorbed 
by another atom in internal state v and located at the 
origin. According to FD statistics two fermions with the 
same quantum numbers repel each other and FD statis- 
tics forces a regular spacing between the atoms. The 
optical interactions are dominantly generated at small 
interatomic distances and the corrections to the suscep- 
tibility due to the second term in Eq. ( |36| ) correspond 
to inhibitited multiple scattering of light resulting in sup- 
pressed diffusive radiation. In the absence of a superfluid 
state FD gas exhibits a dramatic narrowing of the absorp- 
tion linewidth for coherently scattered light jl||,[l6| . 

The first term in Eq. ( pq ) represents the quantum- 
statistical corrections to the reabsorption process be- 
tween atoms in different internal states due to the two- 
particle coherence. This term is nonvanishing only in the 
presence of a superfluid state. Because the total spin of 
an interacting atom pair in Eq. ([!]) is an integer, the pairs 
behave as bosons According to the Bose-Einstein 
statistics two bosons attract each other and the BCS 
pairing favors small interatomic spacing. This results in 
enhanced optical interactions and incoherent scattering 
of light. 

The electric susceptibility exhibits a Lorentzian line 
shape. The optical line shift S and linewidth T for the 
atomic sample are obtained from Eq. (J35|) 



r/ 7 = i 



Gtt 



Ini 



(39) 



(40) 



The collisional line shift <S co i, which results from S2 
[Eq. (p7f)1, is generated by the s-wave interactions. It 
depends on the BCS order parameter A: 



S, 



■col 



p{Ug 



,P2(r ],v i) 



( A 



\TlUgP 



(41) 



The first term in the optical lineshift ( p9| ) corresponds 
to the Lorentz-Lorenz local-field correction. As far as the 



s-wave interactions can be considered local on the scale 
of the optical wavelength in Eqs. ([j]) and ((2^) also the 
line shift 6> co i may be considered as a local-field shift. In 
that case the local-field shift due to the s-wave scattering 
in Eq. ( |39"| ) is larger than the Lorentz-Lorenz shift for 6 Li 
if 



7 £ 210 



1 



( A 



\ flUgP 



[Clg Qgej 3 —1 

— — fim s 
a \ 6 



Here A denotes the wavelength of the incoming light and 
ao is the Bohr radius. Because (A/hu g p) 2 is expected 
to be of the order of one |J, the local-field shift could 
strikingly depend on the BCS order parameter A. The 
collisional line shift was recently observed for a hydrogen 
BE condensate by using a two-photon IS — > 2S spec- 
troscopy EL 

If the the effective range r u of the triplet s-wave po- 
tential in Eqs. ([!]) and fl2C| ) is very short, r u <C 1/fc, the 
resonant dipole-dipole interactions may suppress the ef- 
fect of the s-wave scattering on the line shift just as they 
cancel the effect of the polarization self-energy . How- 
ever, for a metastable state, 7 _1 may be large on the time 
scale of the atomic interactions. In that case the colli- 
sional shift could be observable even for very small r u . 

To calculate the linewidth ( |40"| ) and line shift (^) from 
integral (36) we need to evaluate the spatial correlation 
functions by using Eqs. (|J) and (|4|). For instance, the 
expectation value for the anomalous correlation function 
reads 



(V>i(r)V T (0)> 



V ^ 

k 



2£ k 



(1 



"ak 



)• (42) 



The chemical potential is solved from p v = p v (p). Here 
(O)-0t (0)) = —A/(hu g ) is ultraviolet-divergent, re- 
sulting from the assumption of the contact two-body in- 
teraction in Eq. (Q). This interaction is momentum inde- 
pendent and it is not valid at high energies. To estimate 
the pairing a standard procedure is to remove the high- 
energy divergence by introducing a high-momentum cut- 
off k c . Nevertheless, we find that the optical linewidth 
of a FD gas [Eq. (flof)] is finite even without any high- 
momentum cutoff. This is because the dipole radiation 
already involves a high-frequency cutoff |3l| ] that regu- 
larizes small r behavior. The contribution to the op- 
tical line shift [Eq. (p9|)] from the integral ([36|) is not 
finite. The lower limit of the integral di verg es logarith- 
mically. Although the radiation kernel (15t) involves a 
cutoff |n| , the Lamb shift is not treated rigorously. The 
small-distance singularities of the dipole radiation kernel 
may be regularized by introducing explicit regularization 
factors [|42j. However, for the present purposes we may 
at least obtain an estimate for the shift by using the cut- 
off k c = 1 /r u in the anomalous correlation function d42| ) 
with the realistic value r u = 100ao of the triplet s-wave 
potential [9|. 

In Fig. 2 we have shown (a) the absorption linewidth 
and (b) the line shift for coherently scattered light from 
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Eq. ( p9| ) without the collisional shift, i.e., by assuming 
U g = u ge in Eq. (|3^), for A = 900 nm and for the value of 
the s-wave scattering length of 6 Li, a g — — 2160ao 0. In 
(a) the solid line represents the linewidth in the absence 
of the superffuid state (A = 0). The line narrows as 
a function of the density already at very low densities 
. The presence of the superfluid state broadens the 
optical linewidth (the dashed line). For the BCS state, 
even without the collisional shift, also (b) the optical line 
shift is increased. 

It is interesting to emphasize that the optical linewidth 
is almost independent of the high-momentum behavior 
of the anomalous correlation function ( fl2|). This can be 
seen by introducing a cutoff k c in Eq. ([f2|). We found 
jl7j that the optical linewidth is almost independent of 
the cutoff from k c = oo to k c — l/(500ao) indicating 
that the exact short-range behavior of the two-body s- 
wave potential is not very crucial for the value of the 
linewidth. Furthermore, the contribution of integral j3^ ) 
in the close neighborhood of the origin to the linewidth 
is vanishingly small, and therefore the effects of the BCS 
pair correlations do not result from short-range correla- 
tions. 




FIG. 2. The optical (a) linewidth and (b) the line shift for 
the electric susceptibility, in the absence of the collisional line 
shift, as a function of the atom density per cubic optical wave 
number of the driving light. The scaling of all the variables 
is linear. The solid line represents the optical response in 
the absence of a superfluid state (A = 0). The BCS pairing 
broadens the resonance line and increses the line shift already 
at low atom densities. 



D. Summary remarks 

We calculated the leading quantum-statistical and col- 
lisional corrections to the standard century-old column 
density result for the electric susceptibility. In Eq. ( |35| ) 
the corrections to the susceptibility are encapsulated in 
the two parameters £i and £2. Here £2 represents the 
collisional line shift due to the s-wave interactions and £1 
position correlations between different atoms. The sus- 
ceptibility was obtained by means of the decorrelation ap- 
proximation ([32j) which neglects all the repeated photon 
exchange between the same atoms. These are the micro- 
scopic mechanism for the collective optical linewidths and 
line shifts jn). Therefore, £1 is a direct consequence of 
the quantum-statistical correlations; for an uncorrelated 
atomic sample, with pi(vv, r'er) = p v p a , £1 vanishes. In 
the case of uncorrelated atoms the lowest-order correc- 
tion to the optical linewidth results from the collective 
light scattering p6| , p)|j40[ . This correction is proportional 
to atom density. The quantum statistics is different be- 
cause the correlation length itself depends on the density. 
For instance, in the low-density limit for an ideal FD gas 
we obtain fUSfl : £1 oc (p/k 3 )(£^k) oc p 2 / 3 /k 2 . Hence, at 
least in the absence of a superfluid state, Eq. ([55]) not 
only represents the lowest-order correction to the suscep- 
tibility in terms of microscopic optical interaction pro- 
cesses between the atoms, but it also correctly generates 
the leading low-density correction. 



IV. BRAGG SPECTROSCOPY 

A. Diffraction of atoms 

In this section we consider diffraction of atoms by 
means of light-stimulated transitions of photons between 
two nonparallel laser beams. When an atomic beam 
interacts with a periodic potential formed by a stand- 
ing light wave, it can Bragg diffract, analogous to the 
Bragg diffraction of x rays from a crystal [^3| . The Bragg 
diffraction has been experimentally used as a technique to 
split a BE condensate |^3| and as a spectroscopic method 
to probe the density fluctuations of a BE condensate 

urn- 

In an nth order Bragg scattering process photons are 
absorbed from one beam and stimulated to emit into the 
other n times [^3 24 



Two different momentum states 
are connected by a 2n-photon process. The change of the 
energy of a photon upon the scattering satisfies E = nhuj, 
where uj stands for the frequency difference of the two 
lasers with wave number k. The fractional change of 
frequency upon scattering is assumed to be negligible, 
\uj\ <C kc, so that one finds the familiar relation between 
the change of the wave vector An and the scattering 
angle 9, 



\Ak\ = 2nfcsin(6>/2) 



(43) 







B. Dynamical structure function 

In Bragg spectroscopy, the two intersecting laser beams 
create a moving standing wave with a periodic intensity 
modulation I(rt) = Jcos(Ak, • r — Lot) [ p5| . The in- 
tensity modulation creates an optical potential V(rt) = 
Vq cos(Ak ■ r — Lot) which couples to the local number 
density of ground-state atoms. The dependence of Vq on 
the atomic level scheme and on light polarizations is ana- 
lyzed in Ref. fl . We consider a situation where the inter- 
nal sublevel of the ground-state atom does not change in 
the scattering process. The corresponding Hamiltonian 
density reads 



U B = Vq cos(Ak ■ r - wi^LMVvM 



(44) 



According to Fermi's golden rule the excitation rate is 
then 2ir/h(y /2) 2 S(AK,Lo) pHh where S(Ak , to) is the 
dynamical structure function Ej| 



5(AK, W ) = i^e-^ T |<i|p(AK)|/)| 



xS(huj + Ef- Ei) . 



(45) 



Here Z denotes the grand partition function and the ex- 
pectation value of the density fluctuation operator, 



/%) = £ / 



(46) 



is summed over all possible final states |/), with the en- 
ergy Ef, and thermally averaged over initial states \i), 
with the energy Ef. By using the completeness of |/) 



and p{— q) = jr (q) we may write Eq. (45) as 



S(An,u>) = -1— V / dtd 3 ri d 3 r 2 e wt t 
2ttH J 



iAK-(ri-r2 



^{^(riOj^^rxO)^^*)^^*)) ■ (47) 
We define the static structure function by 

S{An) = h dujS(AK,Lo). (48) 



The dynamical structure function mirrors the velocity 
distribution of atoms and contains qualitative signatures 
of BE and FD statistics || . It displays the modifications 
of the velocity distribution due to the quantum statistics 
including the Fermi inhibition and the Bose enhancement 
of the scattering process. For the case of two BE conden- 
sates it can also exhibit a dramatic dependence of the 
spectrum on the relative phase between the two conden- 
sates |^8|,[30| , and the Bragg diffraction could possibly be 
used as a technique of measuring the relative condensate 
phase. This is because, due to the macroscopic quan- 
tum coherence of the BE condensates, the uncertainty 
in the initial state of the Bragg diffraction may result in 



a destructive or constructive interference of the transi- 
tion amplitudes. The structure function may also provide 
information about the high-energy quasiparticle excita- 
tions |2j| . Here we study the qualitative signatures of a 
superfluid state in the structure function of a FD gas. 



1. Ideal Fermi- Dirac gas 

First, we consider an ideal FD gas studied in Ref. [||. 
We assume a translationally invariant space. In that 
case the correlation function in Eq. (|47]) depends only 
on r = Y2 — Y\. We are interested in the incoherent scat- 
tering processes corresponding to nonforward directions 
with An ^ 0. In the absence of a superfluid state the 
correlation function in Eq. ( |47| ) with v ^ rj represents 
only coherent scattering events. With v = rj we obtain 

(^ !/ (00)V 9 ,(00)^ J? (rt)^ OT (rt)) 

= p 2 + (^(00)^ I) (^))(V(00)^(rt)) . (49) 

For Ak =/= we obtain the dynamical structure function 
from Eq. © 



k,i/ 



• An 



LO + LO R 



(50) 



Here n^ v = [exp {e^/ksT)/ z + denotes the FD oc- 
cupation numbers and z fugacity. We have also defined 
the effective recoil frequency lor by 



H\Ak\ 2 
2m 



(51) 



Expression ([30|) describes a scattering process in which 
an atom in the ground state v with the cm. state k scat- 
ters to the cm. state k — Ak still remaining in the state 
v. The delta function dictates the energy conservation, 
which coincides with the theory for Doppler velocime- 
try of atoms (44| shifted by the effective recoil frequency 

UR %. 

Classical atoms obey Maxwell-Boltzmann statistics 
and their velocities are normally distributed resulting 
in a Gaussian-shaped dynamical structure function [p|. 
Firstly, FD statistics modifies the velocity distribution; 
even an ideal FD gas at T = exhibits a finite width 
in Eq. (^0|). Secondly, the quantum degeneracy affects 
the scattering processes. The product of the occupation 
numbers in Eq. ( |50| ) indicates the Fermi inhibition: The 
scattering events in which an atom recoils to an already 
occupied state are forbidden by the Pauli exclusion prin- 
ciple. 

It is illustrative to describe the Fermi inhibition in mo- 
mentum space |l^] . At T = the fermionic atoms fill the 
Fermi sphere with v = Q(kp— |k|). The scattering sat- 
isfies Eq. p3) . For the first-order Bragg diffraction, with 
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n = 1, all atoms are scattered out of the Fermi sea, if 
|Ak| > 2kp. Moreover, for small scattering angles, 



sin(0/2) < 



A/d 



6ir 2 p 
|A^P 



1/3 



(52) 



the incoherent atomic recoil events are forbidden. When 
the density is increased, at kp > \An\ the scattering is 
at least partially suppressed to all nonforward directions. 

We consider a situation where the densities of the in- 
ternal sublevels are equal p = p^ = pi = N/V. Here N 
denotes the number of atoms. The dynamical structure 
function for an ideal FD gas was displayed in Ref. || . At 
T = 0, and in the degenerate regime with |Ak| < 2kp, 
the result exhibits a characteristic shape of a wedge con- 
sisting of a linear and a quadratic part. The static 
structure function per total number of atoms may be 
evaluated from Eq. (f3(]|). For |A«| > 2kp we obtain 
S(Ak)/(2N) = 1 and for < \ Ak\ < 2k F 



S(Ak) 
2N 



1 



16k F \AK\ 



12\Ak\ 2 - 



JA4 

kp 



(53) 



For free atoms S(Ak)/(2N) = 1 and Eq. (|53J) describes 
the inhibited scattering. 

In the previous discussion we ignored the Rabi oscilla- 
tion dynamics and considered only the transition rates. 
This is a good approximation when only a small frac- 
tion of atoms is scattered, i.e., when the coupling time 
to is much shorter than the oscillation period to -C I/O 
p8| . Here fl denotes the two-photon Rabi frequency. 
To observe the qualitative features in the spectrum the 
coupling time should be at least of the order of the char- 
acteristic time scale in the spectrum, which according to 
Ref. § is (tURep/hy 1 / 2 . 



2. BCS pairing 

In the presence of a superfluid state the correlation 
function in Eq. (E70 at T = reads (for v / rj): 



(f^00)Vv(00)^(ri)VW r *)> 



gv 
= P 2 



|<Vv(oo)VWrt))l 



(54) 



For the case v = r\ we obtain Eq. ( E 9 [) . Analogously to 
Eqs. (H) we need to evaluate Eqs. fl49Q and (jEJ) by means 
of the Bogoliubov quasiparticles [Eq. (||)] . The dynamical 
structure function for incoherent scattering (Aft ^ 0) in 
the quasiparticle vacuum is 



of quasiparticles separated by the wave vector Ak. We 
note that for a superfluid state in the homogeneous space 
there exists a finite energy gap in the excitation spec- 
trum —huj = £k + -Ek-Are > 2|A|. The corresponding 
expression for the static structure function S(Ak) may 
be obtained from Eqs. ( |55| ) and (liF 




FIG. 3. The static structure function per the total number 
of atoms S(Ak) /(2N) as a function of the atom density p in 
units of cm -3 . The change of the atomic recoil wave number 
upon scattering is (a) 10 3 cm -1 , (b) 10 4 cm -1 , and (c) 3x 10 4 
cm - . The solid line represents the diffraction in the absence 
of a superfluid state (A = 0). The BCS pairing dramatically 
increases the incoherent nearforward scattering already at low 
atom densities. 



S(Ak,u) 



A 2 



4^ 

k 



5(hu! + E k + -Ek-A* 



k"k-AK 



1 - 



1 



;k-A K 



E) 



k-Are 



(55) 



wher e th e quasiparticle energies £k are defined in 
Sec. [IB. Expression (Ba) describes creations of pairs 



In Fig. ^ we show the static structure function 
S(Ak)/(2N) for a FD function of the density 

p for three characteristic values of |Ak|. The s-wave 
scattering length a g — — 2160a . The solid line rep- 
resents an ideal FD gas in the absence of a superfluid 
state. The superfluid state dramatically increases the 
structure function (the dashed line) for nearforward scat- 
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tering. The BCS pairing mixes particles and holes near 
the Fermi surface increasing the number of available scat- 
tering channels. This effect is particularly striking for the 
case of small recoil momentum corresponding to nearfor- 
ward scattering. 

We may also consider situations where the internal 
state of atoms is changed in the scattering process. In 
this case the two ground states \g, |) and \g, J.) are 
coupled through a common excited state by the in- 
tersecting laser beams. For instance, the scattering 
rate for the transition \g, J.) — > \g, f) is proportional 
to (ifjg^(00)ipgi(00)ijjg^(rt)ipgi(rt)) and depends on the 
quasiparticle pairing. 



V. CONCLUSIONS 

We studied the interaction of light with a two-species 
atomic supcrfluid gas. Firstly, we considered the propa- 
gation of light and evaluated the quantum-statistical cor- 
rections to the standard column density result for the 
electric susceptibility. Secondly, we analyzed the Bragg 
diffraction of atoms by means of light-stimulated transi- 
tions of photons between two laser beams. The effects 
of BCS pairing may be understood in terms of enhanced 
incoherent scattering processes resulting in the increased 
optical linewidth, line shift, and static structure function. 
These optical properties could possibly signal the pres- 
ence of the superfluid state and determine the value of 
the BCS order parameter in dilute atomic FD gases. 

One particularly promising candidate to undergo the 
BCS transition and to become a superfluid is spin- 
polarized atomic 6 Li. Atoms in two different internal 
levels can interact via s-wave scattering and the 6 Li 
atom has an anomalously large and negative s-wave 
scattering length a ~ — 2160ao- The hyperfme states 
\m s = 1/2, mi = 1) and \m s = 1/2, m, = 0) of 6 Li have 
been predicted to undergo a superfluid transition at 10 -8 
K with a density of 10 12 cm -3 Here m s and m, de- 

note the electron and the nuclear spin components. 

We assumed a translationally invariant system. A FD 
gas in a harmonic trap may be considered locally ho- 
mogeneous ||, provided that the trap length scale I = 
(H/mui) 1 ^ 2 is much larger than the correlation lengths, 
£f! and The spatial confinement introduces an un- 
certainty in the recoil momentum. In the case of Bragg 
spectroscopy, the coherent scattering is negligible, if the 
change of the wave number of the atoms upon scattering 
is larger than the inverse size scale of the atomic sample 
l/l < \Ak\. 
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